Astronomy Sl Astrophysics manuscript no. ms 


©ESO 2013 


March 21, 2013 





Identification of transitional 6\sks in Chamaeleon witli Herschel^ 



A. Riba^nm B. MeriiP, H. Bouj^, C. Alves de OliveiraP, d. R. ArdilsP, E. PugcP, A. KospaP, L. SpezzP, 
N. L.J. CojP, T. PrustP, G. L. PilbratP, Ph. AndrJ^l, L. MatrPI, and R. VavreP 



o 

<N 



1 ESAC-ESA, P.O. Box, 78, 28691 Villanueva de la Canada, Madrid, Spain 

e-mail: aribas@cab.inta-csic.es 
^ Centro de Astrobiologia, INTA-CSIC, P.O. Box - Apdo. de correos 78, Villanueva de la Canada Madrid 28691, Spain 
^ Ingenieria y Servicios Aeroespaciales-ESAC, P.O. Box, 78, 28691 Villanueva de la Canada, Madrid, Spain 
4 Herschel Science Centre, ESAC-ESA, PO. Box, 78, 28691 Villanueva de la Canada, Madrid, Spain 
^ NASA Herschel Science Center, Cahfornia Institute of Technology, 1200 E. Cahfornia Blvd., Pasadena, CA 91 125, USA 
^ Research and Scientific Support Department, ESTEC-ESA, PO Box 299, 2200 AG, Noordwijk, The Netherlands 
^ European Southern Observatory, Karl-Schwarzschild-Strasse 2, 85748, Garching bei Munchen, Germany 
^ Instituut voor Sterrenkunde, KU Leuven, Celestijnenlaan 200D, B-3001, Leuven, Belgium 

^ Laboratoire AIM Paris - Saclay, CEA/DSM - CNRS - Universite Paris Diderot, IRFU, Service d' Astrophysique, Centre d' Etudes 
de Saclay, Orme des Merisiers, 91191 Gif-sur-Yvette, France 
School of Physics, Trinity College Dublin, Dublin 2, Ireland 

Received 19 December 2012; accepted 11 March 2013 



ON 



ABSTRACT 



O 



> 
00 

cn 
o 



Context. Transitional disks are circumstellar disks with inner holes that in some cases are produced by planets and/or substellar 
companions in these systems. For this reason, these disks are extremely important for the study of planetary system formation. 
Aims. The Herschel Space Observatory provides an unique opportunity for studying the outer regions of protoplanetary disks. In this 
work we update previous knowledge on the transitional disks in the Chamaeleon I and II regions with data from the Herschel Gould 
Belt Survey. 

Methods. We propose a new method for transitional disk classification based on the WISE lljjm - PACS 70 jim color, together with 
inspection of the Herschel images. We applied this method to the population of Class II sources in the Chamaeleon region and studied 
the spectral energy distributions of the transitional disks in the sample. We also built the median spectral energy distribution of Class 
II objects in these regions for comparison with transitional disks. 

Results. The proposed method allows a clear separation of the known transitional disks from the Class II sources. We find six 
transitional disks, all previously known, and identify five objects previously thought to be transitional as possibly non-transitional. 
We find higher fluxes at the PACS wavelengths in the sample of transitional disks than those of Class II objects. 
Conclusions. We show the Herschel 70 /um band to be a robust and efficient tool for transitional disk identification. The sensitivity and 
spatial resolution of Herschel reveals a significant contamination level among the previously identified transitional disk candidates for 
the two regions, which calls for a revision of previous samples of transitional disks in other regions. The systematic excess found at the 
PACS bands could be either a result of the mechanism that produces the transitional phase, or an indication of different evolutionary 
paths for transitional disks and Class II sources. 
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1. Introduction 

Protoplanetary disks surrounding young stars are known to 
evolve over timescales of a few million years from a more mas- 
sive and optically thick phase (Class II objects) to optically thin 
debris disk systems (Class III sources; see ,Williams & Cieza 



predicted by disk evolutionary models (|Hartmann et al. 


1998 


ICalvet et al. 


2005 


Fedele et al. 2010l|Sicilia-Aguilar et al. 


2010 


Spezzi et al. 


2012 


). Another important evidence is found in deep 



20111 



for a recent review on the evolution of protoplanetary 
disks). There are several indications of this evolution with time. 
Infrared (IR) observations of star-forming regions show a sys- 
tematic decrease of the IR flux with stellar age (Haisch et al. 
2QQH IGutermuth et al.|[2004l |Sicilia-AguilS et al.,,2006^ |Cur 



mid-IR spectroscopic observations of young stars with disks that 
show dust grain growth, crystallization, and settling to the disk 
mid-plane. These phenomena are found to be correlated with the 
evolution of the disk structure across two or ders of magnitude in 
stellar m ass ( Meeus et al. 2001 ; van Boeke l et al.|2005||Kessler-| 



rie & kenyon|2009| ). In the optical and ultraviolet, observations 



Silacci et al.,2005,,Apai et al.,2005,,Qlofsson et al.|2009| ). 



show that the disk mass accretion rate decreases with time as 



Send offprint requests to: A. Ribas 

* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



Most of the evolution of protoplanetary disks is driven by 
gravitational interaction and viscosity eff'ects in the disk (^ Pringle 
|1981| ). However, some circumstellar disks show evidence of 
a diff'erent evolutionary phase: they are known as transitional 
disks. Compared to Class II disks, they display a clear dip in 
their spectral energy distribution (SED) at short-mid IR (typi- 
cally around 8-12 yum) and a rising SED with flux excesses simi- 
lar to that of Class II sources at longer wavelengths ( [Strom et al.| 
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T989| ICalvet et al.|[2Q02l [Espaillat et aT]|2Q07l [Andrews et aL 



2011). The dips in the SEDs are usually explained in terms of 
dust-depleted regions and/or cavities in the disks, of typical sizes 
of some tens of AU (see |Merm et al. |201Q( [Andrews et al.|2Qll| 
and references therein). 

Several processes have been proposed to explain these gaps 
and holes: grav itational int eraction wit h a low-mass companion 
( [Bryden et aLlp799i Rice et al [[2Q03[ [Papaloizou et al. 2007 ), 



final maps were created using Scanamorphos ( [Roussel[|2012| ) 
for PACS, and the destriper algorithm included in HIPE for 
SPIRE. These two algorithms are optimized for regions such as 
Chamaeleon, which have bright extended emission. Figure [T] 
shows a three-color composite image of Cha I (blue: 70yum, 
green: 1 60 //m, and red: 250yum). 



photo-evaporation ( Clarke et al.|2Q01[ [Alexander et al. 2QQ6a|b|, ^^^P^^ ^^^^^^'^^ 



or gr ain growth ([DuUemond & Domimk 2005} Tanaka et al. 
[20^1 [Bimstiel et al.[[2012| ). Observational evidence of stellar 
or substellar companions has been obtained in some cases (i.e., 
CoKu Tau4 or T Cha, see jlreland & Kraus|20Q8t[Huelamo eTaLl 



Luhman et"aL] (|2008) and 'Lu hman & Muench[ ( [20Q8| ) presented 



2Q11| respectively). If we were able to distinguish between these 



different explanations would better understand the mechanisms 
that produce the gaps in transitional disks, and the planetary for- 
mation scenario. For this reason, any hint on which process gov- 
erns the transition phase is relevant. 

In this paper, we investigate the contri bution of the far-IR 
data from the Herschel Space Observatory ( [Pilbratt et al.|2QlQ) 
to our understanding of transitional disks. We present a new 
method for transitional disk identification and apply it to the 
sample of Class II objects in the Chamaeleon (Cha) I and II re- 
gions. Section [2] describes the data reduction process, the sam- 
ple selection, and the photometry extraction. Section |3] explains 
the proposed method used in the paper to identify and reclassify 
transitional disks. A more detailed discussion of the sample of 
transitional disks is given in Sect.|4] Section [5] summarizes our 
results. 



2. Observations and sample 

2.1. Observations 

The Cha I and II regions were observed by the Herschel Space 
Observatory in the context of the Gould Belt Survey Andre et al.j 
2Q10| ). These regions are part of the Chamaeleon molecular 
cloud complex that also includes the Cha III cloud. The com- 
plex is located at 150-180pc ( W hittet et al.|1997| ) and is one of 
the most often studied low-mass star-forming regions because of 
its proximity. Cha I has an estimated age of ~2 Myr and a pop- 
ulation of -200 young stellar objects ( [Luhman et al.|2008{|Win-| 
ston et al.|2012 |. Cha II harbors a smaller population (~ 60) of 
young sources ( Young et al.|2005t[Spezzi et al.|2008| ). Because 
of their age and location, these regions are perfect scenarios for 
transitional disk search and study. 

Two sets of observations were used for each region: a first 
set taken in parallel mode, using the PACS (70 and 160,yum 
Poglitsch et ari|20TQ| ) and SPIRE (250 350, and 500yum, Grifj 
fin et al.|2010 ) instruments at a speed of 60'7s, and the 100/im 
PACS band at 20'7s from a second set in scan mode. The 
observing strategy is described in more detail in [Andre et al 
( [2010| ). The total observing time in parallel mode for Cha I 
was ~ 8 hours and 6 hours for Cha II, covering a total area of 
~9deg2 (-5.5 and 3.5 deg^). The PACS 100 yum images cov- 



data sets). Obsids for Cha I are 1342213178, 1342213179 (par- 
allel mode) and 1342224782,1342224783 (scan mode), and ob- 
sids 1342213180,1342213181 (parallel mode), and 1342212708, 
1342212709 (scan mode) for Cha II. 

The data were pre-processed using the Herschel interac- 
tive processing environment (HIPE, Ott 2010| ) version 9. The 



the largest census of young stellar obje cts (YSOs) member s of 
Cha I including Spitzer photometry, and [Alcala et al.[ ( [2008| ) and 
Spezzi et"aL] (|2008^ did the same for Cha II. We selected from 



these studies the sources classified as Class II with known ex- 
tinction values. Since we aim to classify transitional disks, we 
also included T25, flagged as Class III in L uhman et al. (2008 ) 
but later found to be a transitional source in jKim et al.|p009| ). 
We rejected objects with signal-to-noise ratio (S/N) < 5 in any of 
the 2MASS bands to ensure a good photometry estimation and 
coordinates measurement. The final sample of Class II objects is 
comprised of 119 sources. 

To our knowledge, 12 sources in the sample are classified as 
transitional disk candidates in the literature: SZ Cha, CS Cha , 
T25, T35, T54, T56, and CHXR 22E from [Kim et aL] ( |2QQ9l ), 
C7-1 from Damjan ov et aT[(|2007l), CR Cha, WW Cha, and ISO- 



Chal 52 from EspaiTIat et al.|p011 
[Alcala etal.[pQ08| ). 



2.3. Photometry 



and ISO-Chall 29 from 



We extracted Herschel photometry of the Class II sample fol- 
lowing these steps: 



1 . We used the Sussextractor algorithm ([ Savage & 01iver|2QQ7 



ered 2.6 deg^ in Cha I and 2 deg^ in Cha II, and add up to a total 
time of 8 hours and 6 hours, respectively (see [Winsto n et al.| 
[2012[ and Spezzi et al. in prep, for a detailed description of the 



in HIPE to detect sources with an S/N > 5 in the PACS im- 
ages. We then visually checked that no obvious source was 
missing. 

2. We cross-matched the initial sample with the detections in 
the PACS images, using a search radius of 5''. This radius 
was chosen based on the size of the point spread functions 
(PSFs) at these wavelengths (~5.8''x 12.1'', 6.7''x 7.3'', and 
11.4"x 13.4", for the 70, 100, and 160yum bands at the cor- 
responding observing speeds). We note that the background 
emission becomes more sifnificant for longer wavelengths, 
producing false detections because of bright filaments and 
ridges. To avoid possible mismatches, we considered as Her- 
^c/z^/-detected sources only those with counterparts at any 
PACS band. For SPIRE, we found the Sussextractor output 
to be highly contaminated with false detections. Therefore, 
we visually inspected the positions of the detected sources 
individually for these bands. 

3. We performed aperture photometry centered on the 2MASS 
coordinates of each detection. We used the recommended 
aperture radii and background estimation annulus for each 
band (see the PACS point-source flux calibration technical 
note from April 2011, and Sect. 5.7.1.2 of the SPIRE data 
reduction guide). The values for the apertures, inner and 
outer annulus radii (in this order) are 12", 20", 30" for 70 and 
100 yum, 22", 30", 40"for 160 yum, 22", 60", 90"for 250yum, 
30", 60", 90"for 350//m, and 42", 60", 90" for 500 yum. 

4. Since aperture photometry was used, objects close to bright 
filaments or cores are likely to suff'er from contamination. 
Also, given the size of the PSF, no photometric measure- 
ments can be performed for close objects (separation less 
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Fig. 1. Left: Three-color composite image of the Cha I region (blue: PACS70yL/m, green: PACS 160yL/m, red: SPIR E 25 Oyt/m). Circles mark 
the position of transitional disks detected in the Herschel images and classified with the proposed method (see Sec |3.3| l. Triangles show sources 
not fulfilling our selection criteria, and squares represent non-detected sources. Right: thumbnails of the IQfim Herschel maps (50 " X 50 The 
color scale ranges from the median value (background level) to 5 cr over this value (black). In both figures, north is up, and east is left. Note that 
WW Cha is located on a bright core. Article number, page 3 of[TT] 



than ~ lO'O- Therefore we rejected ten detections that 
showed obvious problems in their photometry or enclosed 
more than one object. 

After excluding the transitional disks, the final result of this 
process is 41 Class II sources detected at any PACS band (26 and 
15 in Cha I and Cha II, respectively), nine of them detected also 
with SPIRE. 

We checked that the obtained photometry was consistent 
with that from other map-making algorithms (such as photPro- 
ject for PACS), and found no significant deviation. 

We visually inspected the position of non-detected transi- 
tional disks in the original sample and found that C7-1, CHXR 
22E, and ISO-Chall 29 are not detected at any of the Herschel 
bands. Additionally, ISO-Chal 52 is not detected by Suss extrac- 
tor Sit 70yum, but it is at 100 yum. The object is visually found in 
the 70yum image with a flux of 150 mJy over a background root 
mean square (RMS) of 40 mJy. We therefore included the 70 jim 
flux in our analysis. Source T25 is not in the field of view of 
the lOOyum map, which is smaller than the parallel mode obser- 
vations. Coordinates and stellar parameters for the transitional 
disks in this study can be found in Table [T] 

2.4. Photometric uncertainties and upper limits 

The absolute calibration errors for PACS and SPIRE are 5 % 
and 7 % , respectively (see PACS and SPIRE observer manu- 
als). To ensure a conservative error estimation, we used a 15 % 
error value for PACS and 20 % for SPIRE, taking into account 
that the background emission becomes increasingly stronger at 
longer wavelengths. 

When no source was detected, we computed an upper limit 
calculating the RMS of 100 apertures taken around the source, 
using the same aperture radii and correction factors as for the 
detections. The extracted PACS and SPIRE fluxes for the 12 
transitional candidates in the considered sample are reported in 
Tabled 

3. Identification of transitional disks 

3.1. Photometric selection 

Several selection criteria have been used in t he past to sep arate 
transitional disks from Class II sources. Fang et al.| ( |2009| ) used 
a color-color diagram based on the Ks band and on the [5.8], 
[8.0] and [24] Spitzer bands. [Muzerolle et al.| ( [2QT0l ) proposed 
a classification criterion based on the slope of the SED between 
3.6 and 4.8 yum and between 8 to 24 yum. Cieza et al. (2010) also 
used a color-color diagram, based on the Spitzer photometry at 
3.6, 4.5 and 24 jim. However, all these methods were found to 
suffer from different contamination levels, as explained in |Merm| 
[etaL] ( [2QT0l ). 

There is a high diversity in the morphology of transitional 
disks, hence there are various definitions. However, most of 
them share two common characteristics: (1) they have low or 
no excess with respect to the photosphere up to the At^m-oS or 
the pivot point, usually found around ~ 8-lOyum, and (2) they 
have strong excesses for longer wavelengths (see section 7.1 in 
[Williams & Cieza||2011| and references therein). This is trans- 
lated into a decreasing slope of the SED up to Aum-off. and an 
increasing one for longer wavelengths. 

To identify transitional disks using Herschel photometry, we 
computed two spectral indexes (a)\ one between the Ks band and 
12 yum (aKs-u), and the other between 12 yum and 70 yum (^12-70)- 



log(^iF^^)-log(^2Fi^) 

I0g(.il)-l0g(.i2) ' 



cm ^ s ^ Therefore, 



The spectral index is defined as ' 

where A is measured in jim and in erg s"^ ( 
a traces the slope of the SED in the considered range (a > ^ 
rising SED, a < ^ decreasing SED). This s pectral index has 
been intensively used since its introduction by Lada & Wilking| 
(1984 ) to classify protostars and young objects. 

Figure [2] demonstrates that these two slopes together effi- 
ciently separate the two populations. The separation is clearer 
in the 12-70 yum axis, where au-jo < corresponds to typical 
Class II sources, and 0^12-70 > is indicative of the transitional 
nature of the objects. This separation in the slope between 12yum 
and 70 jjm is an expected feature: for short-mid IR wavelengths, 
the slope depends strongly on the presence of weak excess, or 
on the spectral type of the star if there is no excess. On the other 
hand, the definition of transitional disks itself guarantees a pos- 
itive slope for longer wavelengths. This separation also reveals 
the usefulness of the Herschel data for this classification. As 
a result of the selection method, two disks reported in EspaiT] 
|lat et al.| ( [20TT] ), WW Cha and CR Cha, are not separated from 
Class II objects and we confirm that they do not deviate signif- 
icantly from the median SED of the Class II sources in Cha I 
and II (Fig. |5]). Based on this evidence, we consider them as 
non-transitional. The rest of the transitional disks are properly 
separated from Class II sources. The computed upper limits also 
allow us to classify C7-1 and ISO-Chall 29 as non-transitional 
using this method. 

Interestingly, one Class II source shows 0^12-70 > in the for- 
mer diagram. The object, called ESO-Ha 559, has been recently 
identified as a probable edge-on disk in Ro bberto et al.| ( |2012] ) by 
modeling its SED. Its underluminosity with respect to its spec- 
tral type also supports this scenario. We find this type of object 
to be a source of contamination for this method: edge-on disks 
can mimic the SED of transitional sources. Their geometry will 
cause a high circumstellar extinction level, blocking the light 
from the central star at short wavelengths ([Stapelfeldt & Mon 



eti 1999; Wood et al.|2002l|Duchene et al.|2010t|Huelamo et aL 



|2010| ). The disk becomes optically thin for longer wavelengths 
(> 24 yum) and the emission from the star can pass through the 
disk, resulting in an increase of the flux and hence a positive 
slope of the SED. When their spectral type is known, edge-on 
disks can be identified in Hertzsprung-Russell diagrams, as they 
are often underluminous. 

We also note that this method is not suitable for detecting a 
small subsample of transitional disks called anemic ( Lada et aT] 
|2006), homo logously depleted (|Currie & Kenyon|2009 ), or weak 
excess disks ( [Muzerolle et al. pOlOi They are defined as objects 
with low excess at all infrared wavelengths and show Ofexcess < 0- 
For this reason, they cannot be found with the criterion proposed 
in this work. On the other hand, the rest of transitional disks 
should display 0^12-70 > and hence can be properly separated. 



3.2. Morphological classification 

We checked whether any of the remaining seven transitional 
disks were spatially resolved in the Herschel images. Extended 
emission could indeed indicate contamination by a coincident 
background source, a close by object, or the e xtended back- 
ground emission, as shown in|Matra et al.|p012|). 



Given the estimated distance of 160 pc to the Cha I molecular 
cloud ( jWhittet et al.|[TW7l[Luhman & Muenchl[2008l ), the fuU- 
width at half-maximum of the PSFs for the three PACS bands 
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Table 1. Coordinates and stellar parameters of the 12 transitional disk candidates analyzed in this work. 



Name 


I? A 

K.A./2000 


i^ec./2ooo 


A 

Ay 


opl 


T 




A/r 


D 


Kers. 








(mag) 


Type 






(Mo) 


(J<o) 




CR Cha 


10:59:06.97 


-77:01:40.3 


1.5 


K2 


4900 


3.5 


1.9 


2.6 


1,2,3,4,5,6 


CS Cha 


11:02:24.91 


-77:33:35.7 


0.25 


K6 


4205 


1.5 


0.9 


2.3 


1,2,3,5,6,7,8,9 


SZ Cha 


10:58:16.77 


-77:17:17.1 


1.90 


KO 


5250 


1.9 


1.4 


1.7 


1,2,3,4,5,6,7,8,9 


WW Cha 


11:10:00.11 


-76:34:57.9 


4.8 


K5 


4350 


6.5 


1.2 


4.5 


1,2,3,4,5,6,7,8 


T25 


11:07:19.15 


-76:03:04.9 


0.78 


M3 


3470 


0.3 


0.3 


1.5 


1,2,3,4,5,6,7,9 


T35 


11:08:39.05 


-77:16:04.2 


3.5 


MO 


3850 


0.4 


0.6 


0.5 


1,2,3,4,5,6,7,9 


T54 


11:12:42.69 


-77:22:23.1 


1.78 


G8 


5520 


4.1 


2.4 


1.5 


1,2,3,4,5,6,7,9,10,11 


T56 


11:17:37.01 


-77:04:38.1 


0.23 


M0.5 


3720 


0.4 


0.5 


1.6 


1,2,3,4,5,6,7,9 


ISO-Chal 52 


11:04:42.58 


-77:41:57.1 


1.3 


M4 


3370 


0.1 


0.3 


1.0 


2,3,4,5,7 


C7-1 


11:09:42.60 


-77:25:57.9 


5.0 


M5 


3125 








3,4,5,7,12 


CHXR 22E 


11:07:13.30 


-77:43:49.9 


4.79 


M3.5 


3400 


0.2 






3,5,7,9 


ISO-Chall 29 


12:59:10.14 


-77:12:13.9 


5.57 


MO 


3850 


0.65 




1.85 


13,14 



Re ferences. (l)|Gauvin & Strom] ([1992}; (2)|Espaillat et al.|pOTT]); (3)l Luhman (2007); (4) Luhman & Muench|([2008l); (5)|Manoj et al.|([20TT 
(6)|Henning et al.|(|1993|); (7)ILuhman etaLl(|2008|); (8)|Belloche eraLl(pOrii); (9),Kim et aL| ( ,2009, ); (10) |Lafreniere et al.| ( [2008| ); (ll) |Preibischl 



( |1997| ); (12) |Damjanov et alTl|2007| ); (13) |Spezzi et al.| ( |2008| ); (14) [Alcala et al.| ( [2008l ) 



Table 2. Herschel photometry of the 12 transitional disks in the sample. 



Name 




FlOO//m 


Fl60//m 


F250//m 


F350//m 


FsOO/zm 




(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


Detected sources 


CR Cha 


1.61+0.24 


2.19 + 0.33 


2.74 + 0.41 


2.37 + 0.47 


1.69 + 0.34 


1.09 + 0.22 


CS Cha* 


3.08 + 0.46 


2.82 + 0.42 


2.32 + 0.35 


0.88 + 0.18 


0.38 + 0.08 


0.13 + 0.03 


SZ Cha 


3.88 + 0.58 


3.63 + 0.54 


3.86 + 0.58 


2.85 + 0.57 


1.94 + 0.39 


1.14 + 0.23 


WW Cha 


25.91 + 3.88 


32.32 + 4.85 


27.3 + 4.10 


24.92 + 4.99 


12.44 + 2.49 


6.79 + 1.36 


T25 


0.52 + 0.08 




0.50 + 0.08 


0.20 + 0.04 


0.11+0.02 


<0.10 


T35 


0.38 + 0.06 


0.36 + 0.06 


0.200 + 0.03 


<1.69 


<2.10 


<2.06 


T54 


0.60 + 0.09 


0.77 + 0.12 


0.98 + 0.15 


0.46 + 0.09 


<1.04 


<1.18 


T56 


0.68 + 0.10 


0.57 + 0.09 


0.30 + 0.05 


0.30 + 0.05 


0.30 + 0.06 


0.11+0.02 


ISO-Chal 52 


0.15 + 0.02 


0.15 + 0.02 


<1.07 


<1.42 


<2.06 


<2.04 


Undetected sources 


C7-1 


<0.04 


<0.08 


<0.94 


<1.24 


<1.69 


<2.10 


CHXR 22E 


<0.08 


<0.14 


<1.10 


<1.19 


<1.18 


<0.96 


ISO-Chall 29 


<0.04 


<0.07 


<0.85 


<1.41 


<2.65 


<3.00 



SPIRE photometry is very likely contaminated for this source (see appendix). 



(~ 6'^ 7''and 12 '\ see the official PACS PSF documenQ would 
allow us to resolve structures of ~ 900-2000 AU. It is difficult to 
define an outer radius for protoplanetary/transitional disks, but 
typical values range from some tens to ~ 1000 AU in the most 
extreme cases ( [Williams & Cieza|2011| ). Direct imaging of pro- 
plyds and disks in the Trapezium cluster by |Vicente & Alves 
( 2005| ) showed the size distribution to be contained within 50 
and 100 AU. On the other hand, the R^ parameter (defined as 
the radius where the surface density deviates significantly from 
a pow er law and the disk density declines rapidly, see Williams 
& Cieza,^2011^) has typical value s between 15-230 AU (Hughes 
et al.|20Q81 [Andrews et al.|2009l[2QTQ| ). 

This suggests that none of these sources should be resolved 
in the Herschel images. In each of the PACS band, we compared 
the radial profile of the transitional disks with an empirical PSF 
constructed using clean isolated point sources. Of the nine de- 



tected sources, only T54 was found to be extended, as shown in 
Fig. [3] [Matra et al.[ ( |2012 ) showed that the excess beyond 100 yum 
is likely not related to the source, and therefore not originating 
from a circumstellar disk. This interpretation results in a sub- 
stantial decrease in the IR excess coming from T54, making its 
SED not representative of the characteristic inner-hole geometry 
around transitional disks. The case of T54 shows that one needs 
to verify the origin of the IR excess in protoplanetary/transitional 
disks. 



^ http://herscheLesac.esa.int/twiki/pub/Public/PacsCalibrationWeb/ 
bolopsf_20.pdf 



We found no other extended sources in the Herschel images 
and conclude that all the detected transitional disks have far-IR 
excesses related to the sources. Therefore all but one (T54) of 
the transitional disk candidates in the Cha I and II regions are 
confirmed to be point- sources, up to the resolution of the Her- 
schel PACS and SPIRE instruments. 
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Fig. 2. SED slope between the 12 and VOyt/m (au-io) as a function of 
the SED slope between the ^5-band and (aKs-n)- Transitional 

disks from the literature meeting the selection criterion are marked as 
red dots, green downward arrows are those for which only upper lim- 
its could be estimated. Class II o bjects are black d ots. Blue squares 
are pretransitional disks from Espa illat et al.| ( [20TT] ). There is a clear 
separation between Class II and transitional disks due to the different 
shape of their SED. The single black dot with or > is ESO-Hof 559, 
an edge-on disk. This diagram shows the potential for transitional disk 
classification using the 70 jum band. 
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Fig. 3. Average brightness radial profile (black line) for a point source 
(SZ Cha, left) and an extended source (T54, right) compared to the 
observational PSF radial profile (red line). The error bars are the RMS 
of the values. All considered sources present the same behavior as SZ 
Cha except for T54, whose observed radial profile is clearly above the 
PSF profile, indicating that this source is extended and all others are 
point-like in Herschel images. 



3.3. Transitional disk classification 

Thanks to the new Herschel photometric data, we are able to 
reclassify the already known transitional disks in the Cha I and 
II regions based on the shape of their SEDs. CS Cha, SZ Cha, 
T25, T35, T56, and ISO-Chal 52 show a typical transitional disk 

SED. 

Two objects from |Espaillat et al.| pOllj ) do not fulfill our 
selection criterion, which is tuned to identify clear signatures 



of inne r holes. These sources were selected by [Espaillat et al. 
( 2Q11| ) based on their silicate feature strengths, and therefore the 
different results obtained in this study are not surprising. CR 
Cha shows weak excess up to 2 yum and a typical Class II SED 
for longer wavelengths. The SED of WW Cha does not display 
any decrease in its IR emission, typical for transitional disks. 
Indeed, the He rschel im ages support one of the scenarios pro- 
posed by E spaillat et al. (2011): WW Cha is located on one of 
the cores in Cha I and presumably accretes at a high rate. Fur- 
thermore, it shows a strong excess along the whole wavelength 
range and therefore cannot be considered as a transitional disk, 
but ir more likely a Class II source. The dust-rich environment 
around WW Cha might also contaminate the Herschel photom- 
etry and account for part of the observed IR excess emission. 
These sources are then probably non-transitional. Moreover, the 
morphological analysis of the candidates shows that T54 is ex- 
tended (Fig. [3]) and hence a misclassified object. 

Conclusions for the non-detected sources are more compli- 
cated to draw, and they should be treated with caution, since a 
non-detection does not directly reject a candidate, but could sim- 
ply be due to a sensitivity bias. The computed upper limits for 
C7-1 and ISO-Chall 29 exclude them as transitional disks ac- 
cording our selection criterion. ISO-Chall 29 is a special case: 
the upper limit of 35 mJy in PACS 70 yum is lower than the de tec- 
tion of 56 90+8.63 mJy in the MIPS 70 yum band indicated in AT] 
|cala et al.] (2008), and these two measurements are inconsistent. 
ISO-Chall 29 shows both strong Li absorption and Ha emis- 
sion, which confirms it as a YSO ( [Spezzi et al.|2Q08| . However, 
it is the only transitional disk in the sample with photospheric 
fluxes up to 24 yum in our and the |Alcala et al.| ( |20Q8| ) sample. 
These authors also found it to have the steepest aexcess- This ob- 
ject is therefore probably misclassified in the Spitzer images, as 
strongly suggested by the non-detection in any of the Herschel 
bands and by the outlier nature of the object if the MIPS detec- 
tion is considered. We therefore reclassify it as a non-transitional 
source. We stress that our method is unable to detect transitional 
disks with weak excesses, and deeper observations should be 
made to confirm or reject the presence of disks and holes in these 
systems. In the case of CHXR 22E it is not possible to extend 
the analysis further without making strong assumptions about 
the disk mass and morphology. For this reason, we exclude it for 
the remainder of the study. 

As mentioned in Sect. |3.1[ it is important to note that the pro- 
posed criterion will only select transitional objects with increas- 
ing slopes between 12 and 70yum. This feature is unlikely to be 
produced by grain growth alone (see [Williams & Cieza||2011 



for a review on the topic). As a result, the proposed classifica- 
tion criterion may be biased toward detecting only transitional 
disks with large inner holes produced by photoevaporation, gap 
opening by unresolved companions, giant planet formation, or 
a combination of these scenarios. Physical interpretation of this 
peculiar SED shapes requires detailed modeling, and there is no 
full consensus yet on which physical phenomena can be safely 
attributed to each type of SED (see e.g. Birnstiel et al.||2Q12| ). 
A more detailed analysis of this topic will help to determine the 
real impact of this selection effect. 

From the initial sample of 12 transitional disk candidates in 
the Cha I and II regions we confirm six objects to be transi- 
tional disks, reject five sources by photometric or morphologi- 
cal criteria, and leave one object unclassified since it is not de- 
tected in the Herschel images. These numbers imply a signif- 
icant (~ 45 %) observed contamination level in the transitional 
disk sample considered in this study. Given the small number 
statistics, it is premature to extend this result to other samples. 
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However, this result calls for a revision of the known transi- 
tional disks: if applicable to the whole sample, this contamina- 
tion level would imply a shorter transitional-phase lifetime and 
hence could shed some light on the mechanisms responsible for 
the evolution of protoplanetary disks. 

4. Transitional disl<s in the sample 

We searched for additional photometric values in the literature 
for each of t he transitional objects in our sample. |Gauvin &] 
[Stromal 992)) reported optical photometry for all sources in Cha I 
except for CHXR 22E, ISO-Chal 52, and C7-1. We have queried 
the VizieR catalog service and retrieved add itional dat a for these 
target s from GALEX (Ma rtin et al.|[2005]), 2M ASS ( [Skrutskie 



et al.|2006|) DENIS ( DENIS Co nsortium |2QQ5|), WI SE ( [Wright 
et al.l|2Q1 0D, and AKARI ( [Murakami et al.||2007 |). To avoid 



possible mismatching, we chose a search radius of 1 ''around 
the 2MASS coordinates. We rejected DENIS photometry for 
T25, T35, and ISO-Chal 52, since it clearly disagreed with the 
rest of photometric data. (Sub)millimeter data at 870 yum and 
1.3 mm were also included from Belloche et al. ( 201 1 ) and Hen- 
ning et aL] ( [1993[ ), respec tively. Spitzer photometric measure- 
ments were included from Damjanov et al. ( 2007[),[ Luhman et al. 
( [2QQ8l ), [Luhman&Mu ench(2008[), and[A ^ 



cala et al. (2008 ). We 



also retrieved the S pitzer IRS spectra from the CASSIS database 
( Lebouteiller et al.,201 1 ). The resulting SEDs for all sources are 
shown in Figs.[4|[5| and [6] Thumbnails of the transitional disks 
as seen in the Herschel 70//m images can be found in Fig. [T] 
We note here that the cross- shaped PSF at 70 yum is produced by 
the parallel mode observations, and does not represent resolved 
objects. 

We compared the Herschel fluxes of the sample of transi- 
tional disks with the Class II sources in the same region. For 
this purpose, we first inspected the SEDs of all Class II sources 
(after removing the transitional disks sample). We identified 
and removed one object (111 11 1083-7641574) previously clas- 
sified as_an edge-on disk ( [Luhman & Muench|[2008[ [Robberto 



et al.[[2012) ). It could not be identified in the slope-slope dia- 
gram since it is only detected at lOOyum. The object is present 
in the Herschel images at 70 yum, but was not detected by 5"^^- 
sextractor with the selected threshold. We built the median SED 
of all Class II objects, extinction corrected and normalized to 
the /-band. We also computed SED of the first (< 25 %) and 
fourth (> 75 %) quartiles. Given the low detection numbers for 
the SPIRE bands, we used the lowest and highest values instead 
of the quartiles at those wavelengths. We included photometry 
from the R, /, 2MASS, Spitzer, and Herschel bands. We con- 
sidered all objects detected in each band (regardless of whether 
they were undetected in the other bands). The obtained values 
are given in Table |3] We also found the median SED not to vary 
significantly when only K or M type stars were considered. 

The comparison between the SEDs of transitional disks and 
the median SED of Class II objects in Cha I and II shows a sys- 
tematic diflTerence in the 70-160yum range. The six transitional 
disks found with the selection criterion used in this study display 
a clear excess over the obtained Class II median SED, and all of 
them are over the fourth quartile level (uncertainties for T35 are 
consistent with a flux value below this level). 

Even though this median SED was built with a relatively 
small statistical sample, this result clearlyshows that transitional 
disks are brighter at 70-160yum than typical Class II sources in 
these regions. 

Similar phenomena were already tentatively described by 



Winston et al. (2012) in a preliminary study of the YSO pop- 



ulation of Cha I and by [Cieza et aL] ( [201 1[ ) for T Cha, but here 
we provide consistent results derived from a much larger sample 
of transitional disks. This excess was not previously found by 
large prog rams using the Spitze r Space Telescope, such as cores 
2 disks (Ev ans et al.[2003|[^009[ ). This can probably be explained 
by the lower sensitivity of Spitzer at 70 yum. 

A bias toward the brightest objects could aflTect these results 
in two diflTerent ways: we might miss the faintest transitional 
disks and the faintest Class II sources. In the first case, the Her- 
schel data are enough to classify eleven out of the twelve previ- 
ously known transitional objects in the sample (less than 10% 
objects missed). This suggests that the proposed method does 
not suflTer from a strong bias eflTect. The existence of an unknown 
population of transitional disks not identified with Spitzer cannot 
be rul ed out (although this possibility is unlikely, see |Merm et al.| 
[2010[ ). However, this would not alter the systematic diflTerence 
found at the 70-160yum range between the detected previously 
known transitional disks and Class II objects in these regions. 
On the other hand, the second scenario (e.g. missing faint Class 
II sources) would imply lower values for the Class II median 
SED in the Herschel range, producing an even stronger diflTer- 
ence between transitional and Class II disks. As a result, the 
Class II median SED computed here should be considered as an 
upper limit. 

If the excess at PACS bands is a common feature of transi- 
tional objects, two explanations can be proposed to explain it: 
(1) transitional disks are not the result of the evolution of Class 
II sources, but a parallel evolutionary path, or (2) the discrep- 
ancy between transitional disks and Class II objects is produced 
during the transitional phase (maybe even by the same mecha- 
nism that causes the transitional disk evolution). In this case, the 
piling-up of mass at some point of the outer disk could produce 
the ste ep slope and the excess found at 70 yum ( Williams & Cieza 
[201 1[ ). With the available data it is not possible to favor any of 
these scenarios, so we leave this question open to future studies. 

A larger and statistically more significant sample of transi- 
tional disks and modeling are required to confirm whether the 
diflTerence found at the PACS bands applies to the whole transi- 
tional disk sample, to identify the real cause (or causes) of the 
excess, and to understand whether transitional disks are indeed 
a later stage of Class II objects or follow a diflTerent evolutionary 
path. 

5. Conclusions 

We presented a new method for identifying transitional disks 
based on the slope between the WISE 12yum and PACS lOjim 
bands. We have applied this method to the whole sample of 
known Class II objects in the Cha I and II star-forming re- 
gions. We were able to separate known transitional sources 
from Class II objects, and reclassified five objects as possibly 
non-transitional. This method could produce false positives due 
to the presence of edge-on disks, and Hertzsprung-Russell dia- 
grams should be used to reject underluminous sources. As a re- 
sult, we found an observed contamination level of ~ 45% among 
previously identified transitional disks in these regions. The size 
of our sample is relatively small, and these results cannot be ap- 
plied to the whole transitional disk sample. However, a revision 
of the transitional disk population in other star- forming regions 
is warranted to determine the real contamination level and to ac- 
count for its eflfects. 

We built the median SED of Class II sources in the regions, 
and found significantly higher PACS fluxes in the transitional 
disks compared to it. This excess could be produced during the 
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Fig. 4. SEDs of the detected transitional disk candidates, confirmed by our classification criterion and updated with the fluxes from Herschel. 
Black dots are the dereddened observed values from the literature, downward black arrows are flux upper limits from the literature. Herschel data 
are represented in red (squares for detections, downward arrows for upper limits). Uncertainties are within symbol sizes. The IRS spectra from 
[Manoj et al. ( 201 1 ) (black solid lines) and the photospheres (dashed black lines from the NextGen models from Allard et al.| 2012| ) are also plotted. 
The median Class II SED (blue solid line) and the first and fourth quartiles (blue area) are shown (see also Sect.|4]and Table [3]). 




Wavelength (/xm) Wavelength (/xm) Wavelength (/xm) 

Fig. 5. SEDs of transitional disk that do not fulfill our classification criteria. CR Cha and WW Cha display clear Class II SED. T54 appears 
extended in the Herschel images. Symbols are as in Fig. [4] 

C7-1 CHXR 22E ISO-Chall 29 
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Wavelength (/xm) Wavelength (/xm) Wavelength (/xm) 

Fig. 6. SED of the transitional disk candidates not detected by Herschel for which only flux upper limits could be estimated. Symbols are the 
same as in Fig.|4] The upper limits at 70 //m for C7-1 and ISO-Chall 29 allow us to classify them as non-transitional with our selection criteria. 
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Table 3. Normalized flux densities of the median SED, upper SED (fourth quartile) and lower SED (first quartile) of the Class II objects in Cha 
I and Cha II, after extinction correction. Transitional disks are not included. The number of stars detected in each band is also indicated. For 
comparison, we also include the median value for the six transitional disks detected in this study, although we stress its very low number statistics. 
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1.22 


1.02 
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0.67 
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7.5e-5 


2.8e-6 


9 


SPIRE 500 


5.2e-6 


9.1e-6 


2.3e-5 


6.3e-7 


7 



transitional phase, or be explained in terms of a different evolu- 
tionary path for transitional disks and Class II sources. 

Euture studies of other star-forming regions observed by the 
Herschel Space Observatory will clarify the contamination level 
of the sample of known transitional objects, and will provide 
stronger evidence for a systematic excess at PACS wavelengths 
in transitional disks with respect to Class II sources. 
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Appendix A: Literature review of the individual 
transitional objects detected with Herschel 

Appendix A. 1: CS Cha 

CS Cha was first studied by |Gauvin & Strom] ( |1992| ), who found 
evidence that it harbors a disk with inner holes. It is known 
to be a spectroscopic binary system, as confirmed by 'Guenther' 
et al.| ( |2Q0'7] ) (period > 7 years, minimum mass of the companion 
-O.IM©), although previously Takami et al. (2003 ) suggested 
this option based on the large gap found in its disk. This pre- 
viously unknown feature is probably the reason for the spectral 



CS Cha into a circumbinary disk. The disk has been modeled 
intensively, initially excluding the eff'ect of the binary system 
( [Espaillat "et al.^ 200'7l |2011|) and ev idence of an inner hole of 
~ 40 AU was found. Espaillat et al.| ( |2011| ) also pointed out the 
need for a diff'erent mass distribution in CS Cha compared to that 
of disks around single stars. A more recent analysis by |Nagel 
|et al.| (!2012) also accounted for the binary eff'ect. To reproduce 
the variations found at the IR wavelengths, the model includes 
the emission from the inner disk structure generated by the dou- 
ble system, with a ring and streams of material falling from the 
ring to the circumstellar disks around the individual stars. An- 
other 2MASS source is found at 5'' away. It is 6 magnitudes 
weaker than CS Cha in the 2MASS J band and undetected in 
the rest of the 2MASS bands. Contamination from this object is 
therefore very unlikely. 

CS Cha is located in front of a bright background. There- 
fore, the SPIRE fluxes are very likely underestimated because 
the background emission was probably overestimated during the 
photometry extraction. 



Appendix A.2: SZ Cha 



This source was cataloged as a KO star by Rydgren] ( |1980|) and 
was first identified as a disk with a possible inner gap by |Gauvin 
& Strom| ( |1992| ). |Luhman| ( 2007|) reviewe d its properties, and it 



was lately confirmed by| Kim et al.| p009| ) as a transitional disk. 
It has sometimes been referred to as a pretransitional disk, given 
th e small e xcess found at 3-lOyum. The first modeling results 
by|Kim et al. (2009) suggested an inner disk radius of ~ 30 AU. 



[Espaillat et ai.j ( 12011) modeled this object in detail, noting flux 
variations from IRS spectra at difl'erent epochs on periods shorter 
than three years. These variations are attributed to changes in the 
height of the optically thick disk wall (from 0.006 to 0.08 AU), 
and they do not modify the lO yum silicate emiss ion feature. SZ 
is known to be a wide binary ( |Vogt et aLl|2012| ). A companion 
is found at ~ 5'' (projected distance of 845 AU), which could be 
causing truncation of the outer disk. The contribution of this 
source to the total measured fluxes in this study is likely to be 
negligible, since it is 4 magnitudes weaker than SZ Cha in the 
2MASS J band. However, the possibility of an increase in its 
FIR measurements cannot be excluded. 



Appendix A.3: T25 



T25 was identified as an M3 star by [Luhman et al. p008[ ) and 
was found to be a transitional disk by [Kim et al.[ ( [2009] r The 
lack of IR excess at wavelengths < 8 fim indicates that the inner 
regions of the disk are well depleted of small dust particles. The 
modeling by Kim et al.| (|2009 ) yields an inner radius of 8 AU for 
the disk. It is the only detected transitional object, together with 
T35, lacking the silicates feature at lOyum, another indication of 
an efficient depletion of small particles in the inner disk region. 
T25 has no known stellar companions ( [Nguyen et al.|2012 ). 



type inconsistency found in the literature (Henize & Mendoza v 
W73llAppenzeller 1977; Rydgren 1980; Appenzeller et al.|1983 



Luhman 2004 ). In this study we used the K6 spectral type found 
by iLuhman, ( ,2007 J . The binary nature makes the disk around 



Appendix A.4: T35 

[Gauvin & Strom[ ( [l992| ) classified this source as an MO sta r, and 
it was later identified as a possible a pretransitional disk by [Kim 
[et al.| ( [2009[ ) because it displays weak excess at short IR wave- 
lengths. In this case, the inner disk radius is located at 15 AU 
( [Espaillat et al.[[20TT] ). As in T25, there is no sign of silicate 
emission. The excess at 70yum is lower than in other cases, but 
does not resemble the typical Class II SED. It has no confirmed 
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known stellar or substellar companions ( [Nguyen et al.||2Q12] ). 
However, recent sparse aperture masking observations of this 
source by Cieza et al. (2013 ) showed and asymmetry in its K- 
band flux. On the basis of modeling, these authors found the 
inner disk radius to be -8.3 AU. They were also unable to distin- 
guish between the close-companion scenario and the asymmetry 
being produced by the starlight scattered ofl' the disk itself. 



sion, which aflTected our photometry and, hence, our conclusions. 
The non-transitional nature of this object is also supported by the 
fact that it would be the only transitional disk in our sample with 
no excess emission at 70yum with respect to the median SED 
Class II disks. 



Appendix A.5: T56 

This source was found to be an MO. 5 start inl Gauvin & StronJ 
( |1992| ). |Kim et al.| ( |2009| ) identified it as a transitional disk with 
a inner disk radius of 1 8 AU. As in the other transitional disks in 
this study, its excess is higher than the expected Class II flux at 
the PACS b ands. It has no known bound companions ([Nguyen] 
let al.|2012| ). 



Appendix A.6: ISO-Chal 52 



ISO-Chal 52 is an M4 star ( |Luhm an "200?). [Espaiflat et"aT 



([2011) proposed it as a transitional disk, finding the source to 
be an extreme case among their sample: based on variations of 
its Spitzer spectrum, models require the inner wall height to in- 
crease by about 400 % (from 0.0006 to 0.003 AU). We also found 
it to be an outlier in the sense that it has the flattest SED between 
12 and 70//m. No bound companions are known for ISO-Chal 
52 ([Nguyen et al.|2012[). 



Appendix A.7: CR Cha 

CR Cha is an M0.5 star (Gauvi n & S5om|[T992] ). [Furlan et al.| 
(|2009) found it to be an outlier in their sample when compar- 
ing the equivalent width of the silicates emission and the spec- 
tral slope between 13 and 31 yum: it was beyond the parameter 
space considered in their study. The explanation given in|Furlan[ 
[et al.| ([2009]) is that this source could be a pretransitional disk. 
For this reason, E spaillat et al.^ (2011) included it in their sam- 
ple of transitional disks. In this study, we found this object to 
be compatible with a Class II object. It is also located among 
other Class II objects using the proposed classification method 
(Fig|2]). Therefore, although we cannot completely rule out the 
possibility that this object is in a pretransitional disk phase given 
its strong silicates emission, it would be in a very early stage of 
the transitional phase. 



Appendix A.8: WW Cha 



This source was first classified as a K5 object by Gauvin & Strom[ 



( [1992[ ). It was included in the analysis of [Espaillat et a l. (20111) 
for the same reason as CR Cha, and modeled as a pretransitional 
disk. Comparison with the median SED of the Class II sources 
shows that WW Cha is well above the median. The SED of WW 
Cha resembles a typical Class II object, probably still embedded 
in the core, as suggested by its high extinction (Ay -4.8 mag) 
and its position in the Herschel maps. The dusty environment 
in which it is located could significantly pollute the photometry 
and, hence, our conclusions about this object. 



Appendix A.9: T54 

T54 is known to be a misclassified transitional disk ( [Matra et aL] 
2012| ), and therefore we excluded it from our analysis. The Her- 
schel images show contamination from close-by extended emis- 
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